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ABSTRACT 

We present a comprehensive set of stellar evolution models for Procyon A in an effort to 
guide future measurements of both traditional stellar parameters and seismic frequencies towards 
constraining the amount of core overshoot in Procyon A and possibly other stars. Current 
observational measurements of Procyon A when combined with traditional stellar modeling only 
place a large upper limit on overshoot of a?ov < 1-1- By carrying out a detailed pulsation 
analysis, we further demonstrate, how p- and g-mode averaged spacings can be used to gain better 
estimates of the core size. For both p- and g-modes, the frequency spacings for models without 
overshoot are clearly separated from the models with overshoot. In addition, measurements of the 
I = averaged small p-mode spacings could be used to establish Procyon A's evolutionary stage. 
For a fixed implementation of overshoot and under favorable circumstances, the g-mode spacings 
can be used to determine the overshoot extent to an accuracy of ±0.05 Hp. However, we stress 
that considerable confusion is added due to the unknown treatment of the overshoot region. This 
ambiguity might be removed by analyzing many different stars. A simple non-local convection 
theory developed by Kuhfufi is implemented in our stellar evolution code and contrasted with the 
traditional approaches. We show that this theory supports a moderate increase of the amount 
of convective overshoot with stellar mass of Aaov — +0.10 between 1.5 M© and 15 Mq. This 
theory places an upper limit on Procyon A's core overshoot extent of ~ 0.4 Hp which matches 
the limit imposed by Roxburgh's integral criterion. 

Subject headings: stars: evolution — stars: individual (Procyon A) — stars: interiors — stars: oscillations 
1. Introduction 

1 Visiting Astronomer, Department of Astronomy, Yale 

University One of the most important deficiencies in stel- 

lar evolution theory is the lack of a physically 



correct treatment of convective motions and con- 
vective transport of energy. The extreme dif- 
ficulty of solving the well-known hydrodynamic 
equations on all different scales is a long stand- 
ing mathematical and computational challenge. 
Due to these unsolved complications of an accu- 
rate treatment, convection in stars has been de- 
scribed by mixing length theory (MLT) (Bohm- 
Vitense 1958), a simple phenomenological method 
for treating convection locally. Within the frame- 
work of MLT Roxburgh (1965) and Saslaw & 
Schwarzschild (1965) were the first who estimated 
the penetration of convective motions into the ra- 
diative regions due to a finite convective velocity 
at the Schwarzschild boundary defined by vanish- 
ing acceleration through buoyancy. The result, a 
negligible amount of overshoot of 10 -2 Hp was crit- 
icized by Shaviv & Salpeter (1973) who found ap- 
preciable overshoot within the framework of a sim- 
plified non-local mixing length model. Most stellar 
modelers since have considered overshoot only in 
a parameterized way by extending the core size by 
some fraction aov of the pressure scale height Hp 
taken at the locus of the Schwarzschild boundary. 

In the present work we implement a physically 
better motivated phenomenological model for con- 
vection including overshoot into stellar evolution 
calculations, a non-local convection theory devel- 
oped by Kuhfufi (1986) and we contrast this ap- 
proach with the traditional ones. In both cases we 
are left with free parameters. In MLT only the 
overshoot parameter aov has to be considered for 
the extension of the convective cores. The situa- 
tion is more complicated for the Kuhfufi (1986) 
convection theory (KCT). Nonetheless, we vary 
only one parameter a t that is the most relevant 
for the amount of overshoot. 

Empirically, isochrone fitting to color-magnitude 
diagrams (CMDs) of open clusters have so far 
given the most quantitative results for the amount 
of core overshoot. The underlying hypothesis in 
these studies is that all stars have been formed 
from the same material at the same epoch. Both 
age and metallicity of the cluster arc not known 
and a range of parameters has to be matched 
to observations. However, most uncertainty is 
introduced through contamination of unresolved 
binary systems and variable stars, the need to 
clean the sample from non-member field stars 
(cf., Kozhurina-Platais et al. 1997), the relatively 



small number of stars in the sample and the sys- 
tematics introduced by transforming the CMD to 
the theoretical HRD diagram. Studies of this kind 
have been performed to test the overshoot param- 
eter, e.g., by Prather & Demarquc (1974), Maeder 
& Mermilliod (1981), Pols et al. (1998) and Dc- 
marque, Sarajedini, & Guo (1994) who obtain for 
their best fit a value of aov = 0.23, which can be 
regarded as a canonical value. 

In a different approach, detached and evolved 
eclipsing binaries have been used to constrain the 
amount of overshoot. Ribas, Jordi, & Gimenez 
(2000) studied eight such objects and found an 
overall best fit of aov = 0.25 ± 0.05 and more 
strikingly a clear trend of increased overshoot with 
increasing mass. 

Recently, Aerts et al. (2003); Dupret et al. 
(2004) have used observed seismic constraints to 
deduct the extent of core overshoot in the [3 Cep 
star HD 129929. Based on fitting three identi- 
fied pulsation frequencies these authors claim an 
overshoot parameter of aov = 0.10 ± 0.05 for this 
~ 9 Mq star. While unambiguously identified sin- 
gle frequencies hold the promise of strongly con- 
straining stellar models, we think that many more 
frequencies than free stellar parameters must be 
used for such an analysis. We show in this pa- 
per, that the frequencies are also influenced by 
the adopted overshoot prescription, which con- 
stitutes an additional uncertainty not previously 
taken into account in stellar modeling. 

Our main objective in this study is to show how 
averaged frequency spacings can be used to con- 
strain core overshoot in carefully selected single 
stars like Procyon A . The advantage of this tech- 
nique is that averaged spacings may not depend 
on individual mode identifications as critically as 
single modes do. Results from this approach pro- 
duce weaker constraints on stellar models but on 
the other hand should yield more solid results. By 
extensively analyzing the pulsation properties of 
many different models for Procyon A we show that 
the oscillation characteristics contain very detailed 
information about overshoot and the stellar inte- 
rior. We outline ways to constrain the overshoot 
parameter that may become practical in the fu- 
ture. 
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2. Stellar Evolution Theory 
2.1. General Input Physics 

We use the Yale Rotating Stellar Evolution 
Code (YREC) (Pinsonneault 1988; Gucnther et al. 
1992) in its non-rotating configuration. The four 
classic stellar structure equations, namely Pois- 
son's, hydrostatic, energy transport and energy 
equations are solved with a Henyey relaxation 
scheme in the interior and a shooting technique 
in the outer envelope and atmosphere. The fitting 
point of the interior model to the envelope is set at 
the fixed mass Mg t = 1CP 12 M*. A nuclear reac- 
tion network is solved separately. The separation 
of the chemical part from the structural part can 
be done safely in the main burning stage of hy- 
drogen burning considered here since the nuclear 
timescale and therefore the change in composition 
is long compared to the thermal and dynamical 
timescales involved. Procyon A has a very thin 
outer convective zone of ~ 1CT 5 M* and in these 
outer layers standard MLT is used. 

In the absence of overshoot, core convection 
is also treated in a standard way by determin- 
ing the regions of convective instability with the 
Schwarzschild criterion. In regions where it indi- 
cates convection, the temperature gradient in the 
energy transport equation is set to the adiabatic 
gradient, thereby assuming fully developed turbu- 
lent convection. The superadiabaticity which can 
be estimated to be typically on the order of 10~ 6 in 
the deep stellar layers is thereby neglected and the 
mixing length parameter present in MLT becomes 
irrelevant. The deviation of 10~ 6 from strict adia- 
baticity is too tiny to have any effect on the stellar 
structure or the pulsation frequencies which arc 
mainly sensitive to the sound speed and the mean 
molecular weight gradient. In addition, chemical 
species are mixed instantaneously within convec- 
tive regions justified by the short mixing timescale 
compared to the nuclear timescale under the con- 
ditions prevalent in this study. 

For the material functions, the OPAL equation 
of state (Rogers, Swenson, & Iglesias 1996) and 
OPAL opacities (Iglesias & Rogers 1996) for tem- 
peratures above logT > 4.12 are used, and in the 
low temperature regime logT < 4.0, we use opaci- 
ties by Alexander & Ferguson (1994). In the tran- 
sition region the low temperature opacity tables 
and the interior opacities are ramp averaged. 



The number of grid points can be crucial for cal- 
culating correct convective core sizes and pulsation 
frequencies. YREC inserts and deletes grid points 
according to specified criteria on structural vari- 
ables (e.g., pressure). With the help of some test 
calculations we find that ~ 900 grid points are suf- 
ficient for both core size and pulsation frequencies. 
Our models generally consist of about 1800 grid 
points, twice as many as minimally required. The 
innermost mass shell is taken at m\ = 10~ 5 M* 
and we verified that this is about an order of mag- 
nitude lower than needed to produce correct pul- 
sation frequencies. 

A note on terminology. Zahn (1991) proposes 
to distinguish between the term overshooting and 
penetration since they have a more precise mean- 
ing in fluid dynamics. In this nomenclature, over- 
shooting is reserved for inefficient penetration, 
that does not alter the stable temperature gradi- 
ent whereas subadiabatic penetration refers to con- 
vective heat transport efficient enough to establish 
a nearly adiabatic temperature gradient. In this 
paper, we prefer to use overshoot as the superordi- 
nate term for all penetration beyond the stability 
boundary regardless of its efficiency. However, we 
refer to penetration that does not alter the tem- 
perature gradient as overmixing. 

2.2. Overshoot Prescriptions 

There exists vast theoretical and experimen- 
tal evidence for overshoot in geophysical fluids, 
laboratory experiments and computer simulations 
(Zahn 1991). With a simple theoretical model uti- 
lizing the usual scaling laws of thermal convec- 
tion Zahn (1991) demonstrates that the subadia- 
batic penetration above a convective core should 
amount to a substantial fraction of the core ra- 
dius in the interior of stars. If correct, overshoot 
would strongly influence the evolution of stars for 
its structure, chemical enrichment, lifetime and 
possibly its fate. 

It has long been recognized that MLT, still com- 
monly used in almost all stellar evolution calcu- 
lations, being a strictly local theory, cannot ac- 
count for the phenomenon of penetrative motions 
beyond the classical boundary of convective cores. 
Many non-local convection theories available from 
geophysics or other fields have been largely ignored 
for two reasons. First, many of them are diffi- 
cult to implement within the numerical scheme 
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adopted in stellar evolution and probably more 
importantly contain short lived timescales of the 
order of the Brunt- Vaisala frequency thus ren- 
dering it impossible to follow stellar evolution on 
timescales comparable to the thermal and nuclear 
timescales. Secondly, many of those theories con- 
tain more than one free parameter giving argu- 
ment to the view that more elaborate descriptions 
may be more correct but nonetheless exhibit ar- 
bitrary degrees of freedom. In this spirit, and 
partly supported by the success of the one pa- 
rameter MLT, all actual stellar evolution calcu- 
lations employing overshoot, extend the boundary 
of the convective core by some fraction aov of 
the local pressure scale height taken at this very 
boundary. The extended region is treated differ- 
ently by different authors, some only taking ad- 
ditional mixing of chemical species into account 
while others also alter the temperature gradient 
in this region. While this approach being simple, 
it cannot be regarded as an adequate method to 
capture the physical phenomenon underlying core 
overshoot and it is not surprising that the value 
of aov is disputed not only over its correct value, 
but also whether the same value should be taken 
at the bottom and the top of convective regions, 
and whether higher values should be used for stars 
with higher mass (Ribas et al. 2000). 

We feel it is time to advance the prescription 
of overshoot employed in stellar evolution calcula- 
tions. We choose to implement a one dimensional 
non-local convection theory developed by Kuhfufi 
(1986) in the simplest possible way. Note that 
KuhfuB (1986) derives an equation for the specific 
turbulent kinetic energy by proper spherically av- 
eraging of the first order perturbed Navier-Stokes 
equations. The principal difficulty that arises in 
all phenomenological theories of convection is to 
model the (unknown) correlation functions of the 
fluctuating quantities. This is done within the 
framework of anelastic and diffusion type approx- 
imations by introducing free parameters for every 
term arising in the equation. Thus, we are left 
with five free parameters. KuhfuB (1986) fixes two 
of them by matching the convective velocity and 
the convective flux with the corresponding MLT 
values. We do not alter those in this study. The 
remaining parameters consist of a mixing length 
parameter omlt (analogous to MLT) and an over- 
shoot parameter a t (analogous to aov)- The third 



free parameter does not appear in the original 
KuhfuB (1986)-theory and was first introduced in 
it by Wuchterl & Feuchtinger (1998), a param- 
eter (3 r that essentially limits the local pressure 
scale height Hp to a more meaningful geometri- 
cal length scale in the center of the star where Hp 
formally goes to infinity. We use the convection 
model by Kuhfufi (1986) for several reasons. His 
model consists of only one equation making it rel- 
atively easy to implement in stellar evolution (Sec- 
tion 2.3), the equation is derived from the proper 
hydrodynamical equations and we trust it to be 
physically meaningful for conditions prevalent in 
stellar cores. In the stationary, strictly local limit 
the cubic equation of MLT is retained when the 
Lcdoux criterion is employed. And finally, the 
KuhfuB (1986)-theory gives the same qualitative 
behavior for the temperature stratification in the 
overshoot region as derived by Zahn (1991). 

2.3. Implementation of Overshoot in YREC 

2.3.1. Original Treatment 

Core overshoot in YREC has been previously 
treated in a way most widely used in stellar evo- 
lution codes. In a first step, the boundary of the 
convective core is determined by the Schwarzschild 
criterion, i.e., regions where V ra d > V a d are la- 
beled as convection zones. Since superadiabatic- 
ity is small in the interior of stars, as discussed 
before, the actual stellar temperature gradient is 
then set to the adiabatic gradient V ac i within the 
convection zone. In a second step, a new bound- 
ary is determined by adding a fraction «om of the 
pressure scale height to the boundary at radius r s 
determined by the Schwarzschild criterion 

r new = r s + a M Hp (r s ) (1) 

where Hp is taken at the Schwarzschild boundary. 
Finally, all chemical elements are mixed homoge- 
neously within the new boundary radius given by 
r ncw . By construction, only the mixed region is 
extended whereas the temperature stratification is 
not affected directly within the overshoot region. 
This approach is most correctly termed overmix- 
ing. We label models using overmixing with CT1 
and call the overshoot parameter in this case «om- 
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2.3.2. Simple Improvement on Previous Treat- 
ment 

It has been pointed out by Zahn (1991) that the 
temperature stratification in the overshoot region 
is rendered almost adiabatic due to efficient heat 
transport in this region and that this allows for 
more overshoot than would be possible with the 
original stratification. The solutions calculated 
from Kuhfufi (1986)'s model also agree with this 
assessment. In order to mimic this effect, we op- 
tionally alter the treatment of overshoot in YREC 
such that the temperature gradient is set to the 
adiabatic gradient also in the overshoot region up 
to r ncw . In the course of evolution, the convective 
core size will differ from the previous approach 
when adopting the same overshoot parameter. It 
is therefore important to make clear which over- 
shoot prescription was being used when discussing 
values of the overshoot parameter. We do so by 
labeling models using this approach with CT2 and 
denote the overshoot parameter with aov- 

2.3.3. Implementation of Kuhfufi (1986) theory 

The basic equation of Kuhfufi (1986)'s model is 
an equation for the time-development of specific 
turbulent kinetic energy, w [ergg -1 ], 
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where D/Di is the Lagrangian time derivative, p 
is density, P t turbulent pressure, r stellar radius, 
Hp pressure scale height, V a d adiabatic gradient, 
1/A = l/(aMLT-Hp) + l/(/3 r r) geometrically lim- 
ited mixing length scale and Eq viscous energy 
dissipation. The first term on the right hand side 
of this equation is the driving term for convection 
and it is proportional to the convective flux, j a , 
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where a s is the turbulent driving parameter, Cp 
the specific heat at constant pressure, and in the 
case of an ideal gas with radiation pressure the di- 
mensional parameters S and (p take on the values 
S = (4 - 3/3)//?, 13 = Pgas/P and ip = 1 respec- 
tively. The turbulent flux j& is essentially con- 
trolled by the buoyancy term and the molecular 



weight gradient in brackets. This term is posi- 
tive (negative) in regions where the Ledoux crite- 
rion would indicate convective instability (stabil- 
ity) giving rise to a source (sink) term in equa- 
tion (2) thereby creating (destroying) turbulent 
kinetic energy. Note that the temperature gra- 
dient V is not known a priori. It is rather a quan- 
tity that must be determined from the convection 
theory (see Appendix A). The second term on 
the right hand side of the basic equation (2) is a 
dissipation term with cd being the dissipation effi- 
ciency parameter and it acts always as a sink. This 
term controls the superadiabaticity in the station- 
ary limit of fully developed convection. 

All terms discussed so far are functions of vari- 
ables defined locally at each mass shell of the star. 
If we were to restrict ourselves to only these vari- 
ables and assume time independence, the convec- 
tive core size would match exactly the size deter- 
mined by the Ledoux criterion or, when neglect- 
ing V M , by the Schwarzschild criterion. Indeed, 
the only deviation from MLT would be a different 
superadiabaticity, but since it is very small in the 
core it does not make a difference for the stellar 
model whether the temperature gradient V equals 
V ad + 1(T 5 or V ad + KT 6 . 

In view of its relevance to the convective core 
size the novel part of Kuhfufi (1986) theory stems 
from the non-local term, the third term in equa- 
tion (2) where jt is given by 
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Since this term is proportional to the second spa- 
tial derivative of the specific turbulent kinetic en- 
ergy it leads to overshoot (increase of convective 
region beyond Schwarzschild boundary) in regions 
of positive curvature and to undershoot (decrease 
of convective region below Schwarzschild bound- 
ary) in regions of negative curvature. In the case 
of core convection only the increase of convective 
core size has been encountered in our studies. The 
free parameters a s and cp, arc factors that deter- 
mine the values of convective velocity and convec- 
tive flux. As Kuhfufi (1987) demonstrates, they 
can be exactly matched to the ones of MLT by 
setting: 
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We use these standard values throughout this 
study but it must be stressed that this choice is 
not obligatory and with the kind of methods de- 
scribed in this paper it will hopefully prove possi- 
ble to derive those parameters from observations 
as demonstrated here for the overshoot parameter 

We neglect the second term on the left hand 
side of the turbulent kinetic energy equation, since 
we are only concerned with the hydrostatic main- 
sequence phase in which the star is thermally re- 
laxed and Dp/Dt is close to zero. For stars that are 
not in thermal equilibrium, e.g., in the pre-main 
sequence phase or in the contraction phase before 
the onset of helium burning, this term must be 
taken into account. We also neglect the last term 
in equation (2), the viscous energy dissipation. 
This term is only important for high velocity flows 
which are not encountered during the star's main- 
sequence lifetime. These simplifications yield 
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We aim at finding a solution for u) in the stationary 
limit D<I>/Di = and we show in Appendix A that 
equation (8) can be cast into the following form: 
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where A-F are functions of the known stellar 
background model. Without the non-local terms, 
finding a stationary solution is simple, similar to 
MLT, one has to solve a cubic algebraic equa- 
tion at every mass shell. However, with the spa- 
tial coupling an analytic expression is not obvi- 
ous. Therefore, we solve the non-linear, time- 
dependent equation (9). There is just one further 
complication. In radiative regions equation (9) 
does not have a stationary solution. One might 
be tempted to set ui = in those regions which 
is a trivial solution to equation (9). However, the 



zero solution branch is completely separated from 
the finite solution path and therefore ui will always 
remain zero. In reality, some finite turbulent ki- 
netic energy is always maintained due to thermal 
noise. Technically, this problem can be overcome 
by adding an artificial term to the right hand side 
of equation (9): 



+ M (10-10)3/2 

to 



(11) 



This artificial term can be motivated based on 
the physical picture that it maintains a noise, or 
seed convection underground of u) = 1CT 10 in ra- 
diative regions. It must be strongly emphasized, 
that seed values about two orders of magnitude 
higher start to falsify the solution. Order of mag- 
nitude lower values are possible in principle but 
one quickly runs into problems related to lack of 
numerical precision. With this carefully chosen 
artificial term, we first initialize u> to some ar- 
bitrary value (typically 10~ 10 ) and subsequently 
use an implicit time integration method to reach 
the equilibrium state. We utilize LIMEX (Linear 
IMplicit Extrapolation Method) (Ehrig & Nowak 
2002) for performing the time-integration. The so- 
lution yields the specific turbulent kinetic energy 
ui at every mass shell. We define shells to be con- 
vective, if: 

x Q < 0.1 (12) 

thereby determining the size of the convective 
core. The boundary of the convective region is 
sharply defined by an extremely steep falloff of w. 
In fact, the transition region can be estimated and 
is on the order of 100 km, orders of magnitudes be- 
low the model resolution. Finally, the temperature 
gradient can be calculated from (see Appendix A): 

V = Vad + Sa)(V r ad-Vad) + (l-a;a>) + H 

(13) 

and it remains close to the adiabatic one over the 
complete convection region. 

We refer to models that were calculated with 
the Kuhfufi convection theory with GTS and the 
associated overshoot parameter with a t . 

2.4. Kuhfufi- versus traditional Overshoot 

A first comparison of core sizes between MLT 
and KCT has been given by Kuhfufi himself in 
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his PhD thesis (KuhfuB 1987). For a 15 Mq with 
c*mlt = 1.5 and a t — 0.25 he deduces a core over- 
shoot parameter equivalent to 0.4 Hp above the 
MLT Schwarzschild boundary. Furthermore, his 
analysis gives large overshoot of 0.3 Hp even for 
the lowest mass stars around 1.1 Mq. 

However, as has been pointed out by Wuchterl 
& Feuchtinger (1998) in the context of the KCT 
when applied to RR-Lyrae stars, the pressure scale 
height Hp must be limited where it exceeds the 
geometrically possible length scale. In the case of 
RR-Lyrae stars, the pressure scale height is limited 
by the distance to the stellar surface. 

A similar situation occurs in the deep interior 
of stars. The pressure scale height goes to infinity 
when approaching the stellar center rendering it 
an inadequate measure of scale length in regions 
close to the center. Again, it is physically plau- 
sible to limit the length scale to the geometrical 
distance of the center yielding a length scale A of: 



A 
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V olmutHp 
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(14) 



introducing the parameter /3 r which we set to 1.0 
throughout this whole paper. In a non-local the- 
ory of convection this limited pressure scale height 
in the central regions influences not only the core 
size of stars with small convective cores but it is 
also important for all stellar core sizes and masses. 

With this modified length scale, we repeat the 
comparison between MLT and KCT for two cases, 
one with at — 0.11 and the other with at — 0.30 
both values in accord with stellar evolution of 
Procyon A (see Section 3.3). We determine the 
amount of overshoot - measured in pressure scale 
heights - needed to reproduce the models employ- 
ing KCT. As can be seen in Figure 1 our models 
show less overshoot compared to KuhfuB's origi- 
nal analysis: with a t — 0.3 we find «ov = 0.17 
for a 1.5 Mq star and aov — 0.29 for a 15 Mq 
star. Since Ribas et al. (2000) claim an increase 
of core overshoot for increasing mass by measure- 
ments of eclipsing binaries we also ask whether 
KCT can naturally explain this behavior. In- 
deed, we can see a moderate increase of overshoot 
with mass. Nonetheless, the increase is not pro- 
nounced enough to explain the large overshoot of 
0.6 claimed for the 11.1 M Q star V380 Cyg. By 
increasing the mixing length parameter «mlt the 
increase of overshoot with mass can be enhanced 



slightly but an overshoot of 0.6 remains out of 
reach. 

It is obvious to think rotation could play an im- 
portant role in increasing the core through rota- 
tional mixing in a star of such high mass. Guinan 
et al. (2000) have carried out a detailed analysis of 
this system reaching the conclusion that the slow 
rotation of this star cannot account for substantial 
additional core size increase. The recently derived 
overshoot of aov = 0.10 ± 0.05 for another high 
mass star with 9M Q (Aerts et al. 2003; Dupret 
et al. 2004) seems to indicate that the large over- 
shoot for V380 Cyg may not be a universal prop- 
erty of high mass stars. 

3. Stellar Models of Procyon A 
3.1. Constraints from Observations 

Procyon A and its companion white dwarf con- 
stitute a visual binary system with a 40 yr period. 
The redetermined mass of Procyon A based on 250 
photographic plates of observations between 1912 
and 1995 combined with modern direct measure- 
ments of the angular separation of the pair uti- 
lizing HST Planetary Camera and ground based 
coronagraph data are given by Girard et al. (2000) 
with a derived mass of 1.497±0.037 Mq. The same 
study also determines the parallax of Procyon A 
to be 0'.'2832 ± 0'.'0015 a value different from the 
Hipparcos result p = 0728593 ± 0'.'00088. We take 
the mean value of both independent measurements 
and adopt the difference to the upper and lower 
bound of these measurements as our la error esti- 
mate thereby yielding p = 0728457 ± 070025. The 
larger adopted uncertainties for the parallax lead 
to larger error bars for Procyon A's radius and lu- 
minosity than otherwise expected if we had used 
the quoted errors from the actual measurements. 

The radius of Procyon A can be retrieved us- 
ing the stellar angular diameter data derived from 
optical interferometry — (5751 ± 0705) x 10~ 3 
(Mozurkewich et al. 1991). Taking the center-to- 
limb variation from detailed model atmospheres 
into account Allende Prieto et al. (2002) correct 
this value slightly: = (5748±0705) x lO" 3 . When 
combined with parallax of p = 07 28457 ± 0'.'0025 
and applying the laws of error propagation to the 
three quantities involved one finds the radius 



2p tan(0 Q /2) 



= 2.070 ± 0.026 R c . 







(15) 
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adopting Oq/2 = 959'.'64±0'.'02 (Chollet & Sinceac 
1999). This result is statistically identical to 
the independent angular diameter measurement of 
Kervella et al. (2004). 

The effective temperature is measured with the 
averaged bolometric flux (Fuhrmann ct al. 1997) 
F BO L = (18.20 ± 0.43) x 10~ 6 erg cm- 2 s" 1 : 

T eff = 7400 x (^^y = 6530±50K (16) 

With an absolute distance of Procyon A from 
earth given by the parallax measurement, d = 
3.0857 x 10 18 cm/p(") we derive the intrinsic lu- 
minosity (Lq = 3.833 x 10 33 ergs" 1 ) 

L = AttcPFbol = (7.016 ± 0.21) Lq (17) 

and we use the value log 10 (L/L Q ) = 0.8461±0.013 
throughout this study. 

3.2. Constructing Tracks consistent with 
Observations 

Since parallax p, angular diameter and bolo- 
metric flux i^BOL are measured with completely 
independent techniques and data sets the derived 
basic stellar parameters luminosity L, radius R 
and effective temperature T e g can be jointly used 
to constrain stellar models for Procyon A . There- 
fore we require all stellar models for Procyon A 
to fall into a three dimensional error box given by 
observational constraints for L, R and T c g. With 
help of the identity L = 4ttR 2 o-T^ s one of the con- 
straints can be projected onto a plane given by 
the remaining two. If we look, for example, at the 
i?-i-plane the error box is given by a hexagonal 
shape. Both ler and 2a limits as inferred from 
Section 3.1 are shown (Figure 2). 

Our model calculations start from the ZAMS 
with the default parameters of mass M = 
1.497 Mq, solar hydrogen content X = 0.70, solar 
metallicity Z = 0.018, mixing length «mlt = 1-7, 
no overshoot and no envelope element diffusion. 
We now vary the mass M = {1.423, 1.435, 1.45, 
[1.46, 1.47, 1.48, 1.485, 1.49, 1.497, 1.509, 1.515, 
1.52, 1.53, 1.534], 1.55, 1.56, 1.571}, hydrogen con- 
tent X = {0.68, [0.69, 0.70, 0.71], 0.72}, metal- 
licity Z = {0.016, [0.017, 0.018, 0.019], 0.020}, 
mixing length parameter {[1.5, 1.6, 1.7]} and con- 
vection description CT = {[1, 2, 3]} with the 



associated overshoot parameters aoM = {[0.0, 
0.2, 0.4, 0.6, 0.8]}, a OV = {[0.2, 0.4, 0.6]} and 
a t = {[0.2, 0.6, 3.0]}. The full set of all com- 
binations required us to calculate 14025 different 
models. Of these, 2938 fell in the observational 2a 
error box and 1541 in the la error box. In order 
to explore the sensitivity of our results to uncer- 
tainties in the input parameters, we also study a 
more restricted parameter set given by all combi- 
nations within the values of the innermost brack- 
ets. The restricted set consists of 3564 models, of 
which 712 qualify within the 2a and 354 in the la 
observational limits. Models with overshoot are 
calculated in the case of three different overshoot 
prescriptions CT1-CT3 by varying the overshoot 
parameter relevant to the adopted theory. 

For each stellar track that falls within the error 
box we pick out five models, two at the inner and 
outer 2a boundaries, two at the inner and outer 
la boundaries and finally one close to the center 
of the error box. In instances in which models do 
not cross the la box only two models are picked 
out. On all picked out models, a full pulsational 
analysis is performed and discussed in Section 4. 

All models within the la error box are found 
to be in the main sequence phase of evolution. 
Chaboyer, Demarque, & Guenther (1999) also 
consider a model in the hydrogen shell burning 
subgiant phase. If we allow the full 2a we also 
cannot exclude shell burning phase models from 
stellar tracks. However, subgiant models are only 
found for initial masses of M < 1.47 Mq. It is 
worth noting, that with overshoot, it becomes in- 
creasingly harder to bring a shell burning model 
in accord with the observational (2a) constraints. 
By calculating two additional sets of models with 
qom = 0.02 and aotvi = 0.05 and carefully tak- 
ing into account all possible combinations within 
the remaining parameters we could verify that tiny 
amounts of overshoot of «om > 0.05 completely 
inhibit the shell burning model as a viable possi- 
bility. 

3.3. Constraints on Core Overshoot 

As can be seen in Figure 2 models with no over- 
shoot arc consistent with the observational con- 
straints. At present, we cannot claim evidence for 
or rule out overshoot on the basis of traditional 
stellar modeling. But we can try to provide up- 
per limits and best fit values on overshoot for the 



8 



different overshoot prescriptions. 

Since the stellar luminosity is — to first order 
— a function of mass, a lower mass value for Pro- 
cyon A enforces more overshoot in an effort to in- 
crease the model luminosity to the observed value. 
Overshoot generally causes larger convective cores, 
thereby transforming more hydrogen into helium 
which in the course of nuclear burning leads to 
a higher mass averaged molecular weight p.. Be- 
ing also proportional to p, 4 the luminosity is sen- 
sitively increased by small amounts of overshoot. 

Best fit values are ascertained by adopting de- 
fault stellar parameters and adding overshoot until 
the tracks coincide with the best fit locus defined 
halfway between the radius and halfway between 
the luminosity uncertainties, the latter taken at 
the best fit radius. The best fit values are aoM = 
0.06, aov = 0.06 and a t = 0.03 (Figure 3). 

Smaller mass and cumlt and larger X, Z shift 
tracks to lower luminosities thereby requiring more 
overshoot in order to bring the tracks back into the 
observational error box. We now simultaneously 
alter all four input parameters such that maximal 
overshoot is needed. 

The first set of upper overshoot limits is derived 
by using a mass of M = 1.46 M©, omlt = 1.6 and 
hydrogen X = 0.71 and metallicity of Z = 0.019. 
These input parameters correspond to the set of 
restricted limits Section 3.2. We denote the limits 
by this procedure weak limits and they are aoM < 
0.65, aov < 0.75. (see Figure 4). 

The second set of upper limits, the strong lim- 
its, are found by repeating this procedure for the 
minimally allowed observational mass of M = 
1.423 M©, q;mlt = 1.6 and the highest values al- 
lowed for in this study for hydrogen and the metal- 
licity, i.e., X = 0.72 and Z = 0.020. Thus we take 
the most extreme values in favor of overshoot. The 
strong limits are aoM < 1.18 and aov < 1.13 (see 
Figure 5). 

An interesting property of the Kuhfufi (1986)- 
thcory is, that the core size cannot be increased 
above a certain limit regardless of the overshoot 
parameter a t . For Procyon, the extent of the over- 
shoot region is not substantially increased above 
a t = 5. This behavior is rooted in the KCT and 
we will discuss this in more detail in the next sec- 
tion (Section 3.4). 



3.4. Extent of Overshoot region with KCT 

It is important to point out that the role of 
the free parameter a t is different compared to the 
free parameters qom or aov- Increasing the latter 
two will always lead to more extended convective 
cores — up until the star is fully convective. On 
the other hand, a t is a factor regulating the rel- 
ative importance of the non-local terms, i.e., the 
third term in equation 8. Increasing the non-local 
coupling does not always lead to more extended 
overshoot. In fact, extremely large values for aov 
lead to smaller overshoot extents. 

The run of the non-local term is shown in Fig- 
ure 6 (left panel) for three different and ex- 
tremely high values for a t . The largest convec- 
tive cores can be attained with a t ~ 100 which 
corresponds to an overshoot distance of ~ OA Hp. 
The model with a t = 1000 produces smaller over- 
shoot regions. The reason for this lies in the run of 
the non-local term itself: above the Schwarzschild 
boundary this term is positive and acts towards 
more extended cores but in the inner regions it is 
negative and acts to decrease the overshoot dis- 
tance. The net effect can be studied in Figure 6 
(right panel) where we show the corresponding 
stationary solution of the turbulent kinetic energy 
High values of a t finally decrease the turbulent ki- 
netic energy in the inner part and thereby the ve- 
locity of the eddies at the Schwarzschild boundary 
such that they lack the energy to travel far. 

In both panels of Figure 6 we show how the 
maximum core size derived from the Kuhfufi the- 
ory compares with the well known integral crite- 
rion derived by Roxburgh (1978). As can be seen 
from the figures, the maximally allowed sizes by 
KCT almost exactly match the core size predicted 
by the integral criterion. 

With this property of the more physically moti- 
vated overshoot prescription it would be possible 
to constrain the uncertainties of mass, hydrogen 
content, metallicity and mixing-length parameter. 
Since the scope of this paper is to find observa- 
tional constraints from traditional stellar modeling 
when combined with seismological data we do not 
explore this possibility in more detail. However, 
when comparing the different overshoot prescrip- 
tions in the next sections we must keep in mind 
that models of Procyon A employing the Kuhfufi- 
thcory inherit convective cores with an overshoot 
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distance not larger than ~ 0.4 Hp. 

4. Pulsation Analysis of Procyon A 
4.1. p-modes 

Ground based attempts to measure the p- 
mode large spacings have so far been inconclusive 
(Brown et al. 1991; Bedford et al. 1995; Mosser 
et al. 1998; Barban et al. 1999; Martic et al. 
1999). Of these works Mosser et al. (1998) derive 
a tentative value for the large spacing of 53 /iHz. 
Barban et al. (1999) and Martic et al. (1999) 
claim to have observed excess power and spacings 
of 55 fiHz which would be compatible with the the- 
oretical predictions (see Figure 7). However, all 
these measurements have been restricted to single 
site measurements producing data contaminated 
by significant amplitudes of daily aliases. 

While p-modc frequencies as a function of radial 
and azimuthal order v{n, I) hold the most accurate 
and detailed information about the run of sound 
speed within a star, one cannot hope to match 
the model to all frequencies by blindly search- 
ing the allowed parameter space. This procedure 
would only work if a stellar model matched the 
star as well as our solar models match the Sun. 
Also, inversion techniques cannot be easily ap- 
plied to stars because we are only able to observe 
low azimuthal order p-modes, hence, cannot con- 
struct detailed enough kernels to unambiguously 
span the surface layers. Nonetheless, we can make 
use of averaged quantities which dispense some 
of the rich information contained in the oscilla- 
tion spectrum but still pose additional and use- 
ful constraints on the stellar model. Such infor- 
mation is available from the averaged first- and 
second-order frequency spacings (cf. , Christensen- 
Dalsgaard 1984, p. 11). 

The frequency spacings for Procyon A have 
been calculated and analyzed extensively before 
(Guenther & Demarquc 1993; Chaboyer et al. 
1999). They noted that mode bumping occurs 
in Procyon A and that the irregular spacing as- 
sociated with mode bumping could be used to pin 
down Procyon's evolutionary phase. Independent 
of the oscillation data, they also noted an incon- 
sistency between the astrometric mass of Procyon 
and the mass obtained from fitting their models 
to Procyon's position in the HRD, with the latter 
turning out to be correct (Girard et al. 2000). 



In this paper we discuss what existing obser- 
vational data can tell us about convective core 
overshoot, and we ask how future observational 
data, including seismic data, can help constrain 
the models further. 

For all our models that fall within the error box 
we calculate the I = 0, 1, 2, 3 p-modes from radial 
order n = up to about 35 using the non-radial 
pulsation code JIG developed by Guenther (1994). 
Mode bumping for Procyon A occurs at low radial 
orders between ~ - 16. Mode bumping, also 
referred to as avoided crossings, occurs in more 
evolved stars when the molecular weight gradi- 
ent left behind by nuclear burning leaves a barrier 
that traps g-modes in the core. When the frequen- 
cies of the g-modes increase, as the star evolves, 
they eventually cross through the frequencies of 
the lowest radial order p-modes. In regions of 
frequency overlap the modes mix, perturbing the 
frequencies. Therefore, mixed modes contain ex- 
tremely interesting information on core overshoot. 
Since we are dealing with a resonance phenomenon 
very accurate stellar models are required. One of 
us (D.B.G) is exploring these possibilities in a sub- 
sequent paper. 

The p-modes with radial order greater than 
about 16 show regular spacings in all our mod- 
els. It is this information that is likely the first to 
be detected and measured accurately. 

4-1.1. Large Spacings 

The first order frequency spacings or large spac- 
ings are defined as 

Av(n,l) = v(n,l)-v(n-l,l) . (18) 

It follows from Tassoul (1980) asymptotic theory 
that high radial order modes are approximately 
equally spaced with 

where R is the radius and c is the sound speed. 
Since the integrand is proportional to the inverse 
of the sound speed the large spacings are most 
sensitive to the outer layers of the star where c 
is smallest. The spacings depend on the radius 
and the mass of the star with 8vjv oc 5r/R* — 
SM/(3M lf ). Can we improve on Procyon A's ra- 
dius estimate by making use of the information 
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provided by the large spacings? Indeed we can, 
as seen in Figure 7. Here we plot the large spac- 
ings averaged over the range n — [3 ... 35] ver- 
sus stellar radius for the order I = 0. Each point 
in the diagram corresponds to a calculated stel- 
lar model with CT2. The spread from a straight 
line shows the full set of model uncertainties due 
to different mass, overshoot, metallicity and mix- 
ing length parameters. We verified that models 
corresponding to the other convection prescrip- 
tions, i.e., from CT1 and CT3 do not widen the 
spread. The spread is tight enough to improve on 
the radius by a factor of two for measured aver- 
aged large spacings < 51 ^Hz and > 54 ^Hz. The 
radius measurement in the range 51-54 fiRz can 
be improved upon only with smaller improvement 
factors and it is impossible to do so in the range 
52-52.5 fiRz. The I = 1,2 modes could provide 
additional independent measurements but it turns 
out that there is more than one region where mod- 
els cluster which leads to a wider spread. 

Since the large spacings are sensitive to the 
outer layers we do not expect them to be sensitive 
to different core overshoot. For different orders of I 
we looked at the convective core size as a function 
of the averaged large spacings. The points scatter 
all over the diagram so no correlation between the 
large spacings and the convective core size exists. 
This fact is indeed fortunate if we want to use the 
large spacings for an independent measure of the 
stellar radius. 

4-1-2. Small Spacings 

The small spacings are defined by 

5u(n,l) = u(n,l) - v(n- 1,1 + 2) . (20) 

Again, asymptotic approximations predict equal 
spacings proportional to an integral with an inte- 
grand of r~ 1 (dc/dr) which increases with depth. 
Thus, the small spacings are most sensitive to the 
interior where the integrand is largest. 

Therefore, we want to know whether the small 
spacings could be used to discriminate between 
different convective core sizes. We show the small 
spacings versus the core size, averaged over radial 
orders n = [16 . . . 34] . If we look at the I = 0, CT2 
models in Figure 8 we see that for frequency spac- 
ings > 3.4 ^Hz only models without overshoot can 
be found. In addition, almost all subgiant mod- 
els which can be identified as the ones exhibiting 



zero convective core sizes occupy exclusively the 
region > AfiRz. This result is robust, because 
it does not depend on the overshoot prescriptions 
since only tiny amounts of overshoot inhibit the 
subgiant models completely Thus, the averaged 
small spacings for I — might provide the neces- 
sary measurement needed to discriminate between 
the two evolutionary stages. 

In comparison to the I = modes the I = 1 fre- 
quency spacings are more influenced by the differ- 
ent amounts of overshoot. Here, overshoot leads 
to larger frequency spacings (Figure 9). Due to 
the spreads in parameter space, inferring the over- 
shoot parameter or convective core size is difficult: 
If, for example, a frequency spacing of 5 /iHz is 
measured, «ov could be anywhere between zero 
and 0.6. The situation is more favorable for a 
measured frequency spacing of 8 /xHz where one 
could infer an «ov ~ 0.4 ± 0.2. As will be shown 
in the next section, the g-modes are better suited 
to yield a core size measurement. 

4.2. g-modes 

Although g-modes have never been detected 
conclusively in the Sun (Appourchaux et al. 2000; 
Wachter et al. 2003), it may be possible to detect 
them in Procyon A. The convection zone in the 
sun's atmosphere constitutes a long attenuating 
tunnel for g-modes (Stix 2002, p. 234) in strong 
contrast to Procyon's very thin convective enve- 
lope of ~ 10~ 5 M*. On this basis we hope that for 
existing and planned asteroseismic observations, 
some effort is made to detect g-modes in Pro- 
cyon A . 

4-2.1. Period Spacings 

For all our models that fall within the error box 
we calculate the I = 1,2,3 g-modes from radial 
order n = -2 to about -30(1 = 1), -70(1 = 2), 
— 80(Z = 3) respectively with the non-radial pulsa- 
tion code JIG. In the following, we use the full set 
of combinations, i.e., full input parameter space 
and 2(7 observational error box. As described in 
Section 3.2 we perform a pulsational analysis on 
five selected models within the error box, if pos- 
sible. For one selected model, we show their pe- 
riods as a function of radial order n (Figures 10). 
There are a few gaps in the range of g-modes com- 
puted, an artifact of the automated search tech- 
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nique used by the pulsation program. This is of 
no consequence here since we are here only inter- 
ested in computing the slope of the curves. We 
use a least-square algorithm to fit linear functions 
to the calculated gravity modes. The steepness 
of these lines gives the period spacings which can 
be explained by asymptotic theory (Section 4.2.2). 
Since asymptotic theory predicts a unique period 
spacing any scatter of calculated frequencies from 
a straight line indicates deviations from this the- 
ory. There is some scatter at very low n where 
asymptotic theory does not hold and also where 
mode bumping is likely to occur. In addition, as 
can be seen in Figure 10 but is present in all calcu- 
lated g-mode spectra, the scatter is larger for lower 
I degree modes especially at high radial order. 
This scatter may well contain information about 
stellar structure beyond the information contained 
in the period spacing but we confine ourselves to 
the latter because the scatter is a quantity much 
more difficult to measure. 

4-2.2. Asymptotic Theory 

According to the asymptotic theory (Tassoul 
1980), the period H n t of low degree I and high 
radial order n (n » I ~ 1) is, to first order, given 
by (Unno et al. 1989, p.83) 



(n + l/2- 1/4) 



IT 



with 
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ln(r b ) 



N(r) dlnr 



ln(r cz ) 
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where r cz is the radius of the core convection zone, 
rb the radius at the bottom of the convective enve- 
lope and N is the Brunt- Vaisala frequency. This 
formula establishes a regular period spacing II; 
which is given by 



n, 
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(23) 



As can be seen by this formula, the higher degree 
I modes do not hold additional information about 
the star's structure. The most important contri- 
butions to the integral in equation (23) come from 
the fully ionized inner regions of the star, hence 



the square of the Brunt- Vaisala frequency, A 2 , is 
well approximated by making use of the gradient 
of the mean molecular weight, i.e., 

"4-3/3, 



N = 



9 2 P 



-(V ad - V) + <^V At 



(24) 



where we denote local gravitational acceleration 
g = (Gm) I r 2 , density p, pressure P, ratio of gas to 
total pressure /3 = P gas /P, adiabatic gradient V a d 
and actual stellar temperature gradient V. The 
Brunt- Vaisala frequency gives the frequency by 
which a bubble of gas oscillates vertically around 
its equilibrium position under gravity. The fre- 
quency is therefore proportional to its restoring 
buoyancy force which is given by the spatial en- 
tropy gradient, or alternatively, by V a d — V. An- 
other contribution to the Brunt- Vaisala frequency 
is given by the second term in brackets of equa- 
tion (24): The molecular weight gradient which 
can increase or decrease the frequency depending 
on the sign of V M . 

4-2.3. g-mode Characteristics in Procyon A 

In the case of the late main-sequence star Pro- 
cyon A the Brunt- Vaisala frequency is largest in 
parts of varying molecular weight gradient which is 
left over from the receding core convection zone of 
main hydrogen burning (Figure 11 ). The molec- 
ular weight gradient for mass shells greater than 
about 0.4 M* remain untouched and the buoyancy 
starts to dominate equation (24) . Although buoy- 
ancy constantly increases toward the outer stellar 
layers the drop-off of density reduces its contri- 
butions to the frequency exhibiting a maximum 
around 0.6 M*. These general characteristics are 
true for all our models irrespective of any over- 
shoot, because this basic structure can be found 
in all late main-sequence stars that have convec- 
tive cores. 

4-2-4- What g-modes reveal about overshoot 

With asymptotic theory (Section 4.2.2), we can 
gain further insight on how overshoot influences 
the g-mode period spacings. In order to show 
the differences in the Brunt- Vaisala frequency in- 
scribed by different amounts of overshoot we show 
Figure 12 where we compare KCT overshoot for 
two different overshoot parameters a t — 0.17 and 
a t — 0.30. We recall, that since overshoot in- 
creases the size of the convective region it taps a 
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larger hydrogen fuel reservoir. Models which fall 
into the same error box and exhibit larger over- 
shoot have not used up hydrogen as much and 
therefore the molecular weight gradient is shal- 
lower albeit increasing the period spacing n . Not 
so obvious at first glance but even more important 
is the inner cutoff radius r cz which is larger for 
models with overshoot, reduces the Brunt- Vaisala 
frequency and therefore also increases the period 
spacing n . 

These effects on the period spacings imply, that 
for given and fixed overshoot prescription and in- 
put parameters (mass, metallicity, etc.) a one-to- 
one relation between period spacings and convec- 
tive core size could be made. On the other hand, 
the situation here is more complicated since we 
must take a range of uncertainties into account. 

Here we present the relation between period 
spacings and convective core size in Figures 13- 
21 for all models of the full input parameter 
set that fall in the 2a observational error box. We 
can see a strong albeit not perfect correlation of 
the period spacing with convective core size. This 
can be used to measure the overshoot parameter 
within a given overshoot prescription provided it 
is possible to derive the period spacing from ob- 
servational data. However, within the current un- 
certainties, this measurement is not extremely ac- 
curate. For adopted overshoot prescription and 
given value of the period spacings the derived over- 
shoot parameter can be determined with an ac- 
curacy of ±0.15 .Hp (CT1), ±0.10 Hp (CT2) and 
±0.05 Hp (CT3). Note that although the over- 
shoot parameter at can only be constrained within 
±1.4 the overshoot extend measured in pressure 
scale heights is most favorably constrained with 
CT3. We will see in the next Section 4.3 what un- 
certainties need to be improved upon in order to 
be able to get a more constrained overshoot esti- 
mate. 

Still, we can deduct valuable information about 
the convective core size, especially in certain pe- 
riod regions. Whereas models without overshoot 
occupy the regions with period spacings of 25 — 
45(Z = 1), 21 - 27(1 = 2) and 12 - 19(7 = 3) 
minutes the models with overshoot have generally 
longer period spacings. Thus we find, that irre- 
spective of the different overshoot prescriptions, 
there exists a period regime for models without 
overshoot that is cleanly separated from the over- 



shoot cases. 

For period spacings greater than 45 (7 = 1), 
21(1 = 2) and 19 (I = 3) we get two different 
types of overlaps. First, for given adopted over- 
shoot prescription and overshoot parameter, there 
is a spread due to the input model and obser- 
vational uncertainties. As seen in Figures 13-21 
these spreads constitute islands which correspond 
to a period regime which overlaps with neighbor- 
ing ones. The less overlap the more accurate it 
is possible to determine the overshoot parameter 
from the period spacings. Note that the overlap 
is stronger for the CT1 models when compared to 
the CT2 ones. The second overlap stems from the 
fact that we do not know the correct overshoot 
theory before hand. For example, if we look at 
the I = 2 plots (Figures 14,17) at a period of 31 
minutes with CT1 we get c«om = 0.35 ± 0.15, on 
the other hand, with CT2 we determine «ov = 
0.3 ±0.10. 

Both types of overlaps could be reduced. The 
first one could be reduced by better observational 
data that is able to narrow down the error box 
and the set of input parameters. The most crucial 
observational parameters for this will be identified 
in the next Section 4.3. The second type of over- 
laps that result from the different overshoot pre- 
scriptions is even harder to overcome: when the 
outlined analysis presented here for Procyon A is 
applied to many different, at least equally good 
constrained stars, it might be possible to do so in 
the future. 

4.3. Sensitivity of Results 

In the last section we demonstrated that, for a 
given convection prescription, uncertainties on our 
models imposed by both the observational error 
box and the possible initial input parameter space 
lead to overlaps in the period spacings thereby 
imposing uncertainties on future measurements of 
the overshoot parameter. For upcoming observa- 
tions on Procyon A it might be enlightening which 
parameters need to be constrained the most in or- 
der to be able to determine the overshoot param- 
eter with more accuracy. 

First we want to determine the effects of nar- 
rowing the observational error box. To do so, we 
only allow those models to qualify as viable mod- 
els that fall into the la observational error box. 
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As seen in Figure 22 the period spread occupied 
by the CT2, I = 2, «ov = 0.2 la error box mod- 
els (triangles) is not smaller compared to the full 
model set (squares). Thus, narrowing the error 
box to la does not help in improving the accu- 
racy of the overshoot parameter. 

Next we assess the possible benefits arising from 
constraining one out of the five input parame- 
ters, i.e., mass (Figure 23), hydrogen fraction (Fig- 
ure 24), metallicity (Figure 25) and mixing length 
parameter (Figure 26) . As can be seen in the Fig- 
ures the most crucial input parameters turn out 
to be both the mass and the metallicity. Measur- 
ing those two more precisely can bring down the 
spread by a factor of three. 

5. Summary and Conclusions 

We have constructed a comprehensive set of 
stellar models for Procyon A with varying model 
assumptions on mass, hydrogen content, metallic- 
ity, mixing length parameter, adopted core over- 
shoot prescription and overshoot parameters in an 
effort to single out current and promising future 
observables that constrain the amount and possi- 
bly the nature of core overshoot in stellar interiors. 

In this study, effects of rotation have not been 
taken into account. Procyon A is a slow rotator 
with a measured rotational velocity of v sin i = 
3.2 kms -1 . Rotation is therefore expected to pro- 
duce multiplcts for higher azimuthal modes than 
the fundamental with splittings on the order of 
0.5 /j,Hz. This has to be kept in mind when com- 
paring the predictions made here with actual ob- 
servations. 

The classical treatment of overshoot in stars — 
extending the core size by some fraction of the 
pressure scale height — is merely a measure of 
overshoot and cannot be regarded as an accu- 
rate description of the underlying physics. For 
this basic reason we implement a physically bet- 
ter motivated non-local convection theory devel- 
oped by Kuhfufi (1986). We neglect a number of 
terms present in this theory, foremost its time- 
dependency and some others expected to be only 
important in evolutionary phases that must be fol- 
lowed on a dynamical timescale. More relevant to 
this study is the choice of the five free parame- 
ters present in the theory. Reasonable choices on 
four can be made in analogy with MLT and we 



only vary one parameter, at, which is the most 
relevant to the amount of overshoot. However, 
we stress that all parameters must be eventually 
calibrated, a task that might prove possible in 
the future if the methods presented in this pa- 
per are applied to many stars alike. Regardless of 
the actual adopted numbers, we do find a mod- 
erate increase of the amount of core overshoot 
with increasing stellar mass up to about 15 M Q . 
In comparison with Kuhfufi 's original analysis we 
find smaller amounts of overshoot for the same 
adopted parameters mainly because we limit the 
pressure scale height in the stellar center. We are 
not certain, whether the implementation by Kuh- 
fufi (1987) omitted the molecular weight gradient 
which would constitute an additional reason why 
he found overall larger convective core sizes. 

We also show, that models calculated with the 
Kuhfufi convection theory contain convective core 
overshoot extents smaller than ~ 0.4 Hp regard- 
less of the overshoot parameter a t adopted. This 
behavior is rooted in the non-local treatment of 
the theory: Since overshooting eddies transport 
internal energy farther out it also lowers the in- 
ternal energy in deeper regions of the core. But 
the deep regions provide the driving necessary for 
overshoot to occur. Thus we have two competing 
effects which lead to a limit for the extent of the 
core size and this limit matches the limit imposed 
by Roxburgh's integral criterion. 

The current observational data of Procyon A 
for its basic stellar parameters of mass, lumi- 
nosity, radius and effective temperature constrain 
the amount of core overshoot in our models only 
marginally. Models without overshoot are in ac- 
cord with the observational error box and possible 
input parameters. Our strong upper limits based 
on traditional stellar modeling for overshoot are 
aoM < 1-18 and aov < 1-13. By restricting the 
input parameters of mass, hydrogen, metallicity 
and mixing length to their la levels we also de- 
termine weaker upper limits, i.e., aoM < 0.65 and 
aov < 0.75. The convective core sizes of all mod- 
els with Kuhfufi convection are much more con- 
strained by the theory itself and not the observa- 
tions. 

Procyon A could be in either the late main- 
sequence or the subgiant phase. However, only 
small amounts of overshoot of aoM > 0.05 exclude 
the subgiant models as a viable possibility. We 
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stress that the subgiant phase might be identified 
by the I = averaged small p-mode spacings. The 
identification of Procyon A as a subgiant would in 
turn yield an extremely precise upper overshoot 
limit of aoM < 0.05. 

Additional information might be extracted 
from averaged pulsation data. In certain fre- 
quency regions, the spread of model uncertainties 
is narrow enough to use the averaged large I = 
spacings of the p-modes to improve on the cur- 
rently known radius of Procyon. Improving the 
radius measurement by a factor of two will not 
lead to better estimates on the convective core 
size, since the largest uncertainties reside within 
the mass and metallicity parameters which we 
were able to show in our sensitivity study on the 
g-modes. Certain frequency regions of the aver- 
aged small spacings are exclusively occupied by 
models without overshoot and in addition, the 
subgiant models are separated in the I = spac- 
ings. The p-mode separations due to overshoot 
are too small to give good estimates on the con- 
vective core size. Extremely large overshoot might 
be detected, e.g., qom ~ 0.6 ± 0.2. 

More direct measures of overshoot are hidden 
in the period spacings of g-modes. In compar- 
ison to the p-modes, the separation of averaged 
period spacings due to different extents of over- 
shoot is more pronounced, especially for the I = 1 
modes. When adopting CT2 the best possible de- 
termination of the core overshoot extent is within 
±0.1 Hp. For CT3 the overshoot extent can be 
determined even more accurate within ±0.05 -Hp. 
Again, most models without overshoot and in ad- 
dition the subgiant models are nicely separated 
from the overshoot models. Since the CT3 mod- 
els have restricted core sizes, this theory could be 
falsified for very large period spacings. 

Apart from these described certain regions, con- 
siderable confusion is added due to the different 
possible overshoot descriptions. At the moment, 
neither can we discriminate between the different 
theories nor can we get precise measurements of 
the overshoot parameter based on the averaged 
spacings alone. To improve on this situation we 
must bring down the model spreads by getting 
more precise values for the traditional stellar pa- 
rameters, the most important in this case are the 
mass and the metallicity. The confusion due to 
different convection theories may be removed by 



similar analysis for a couple of equally well selected 
stars. 
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A. Derivation of Implemented Convection Equations 



First, we derive an expression for the actual temperature gradient V which can be given as a function of 
specific turbulent kinetic energy Q. The total flux l/(47rr 2 ) is given by the sum of convective fluxes j a and 
jt and radiative flux 

l 

3q + jt + F Tad . 



Airr 2 

With the definition for the radiative gradient V ra d and the flux F ra d 

Hp I 



V ra( J 
-^rad 
^-rad 



fc rad T 47rr 2 
16ctT 3 



3np 



we can eliminate the luminosity I and the radiative flux F ra d yielding 

H P 



kradT 



' rad 



V . 



With the expressions for j Q (see eq. [3]) and j t (see eq. [5]) and the following abbreviations: 



we can solve for V yielding: 

V = Vad + %Q (Vrad 



Vui) + (-i -- [G^ 



H 



(Al) 

(A2) 
(A3) 
(A4) 

(A5) 
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&rad 


(A6) 
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a s C P 
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d = 


a t H P f 4Trr 2 p\ 
T V Am J 


(A9) 


F = 


ab 


(A10) 


G = 


d/b 


(All) 


H = 


c/b 


(A12) 




1/(1 + Flu^ 2 ) 


(A13) 



(A14) 



This expression can now be substituted into the original equation for the turbulent kinetic energy (see eq. [8] ) 
so we can arrive at an equation in which the coefficients depend only on known stellar parameters. With the 
additional abbreviations: 



I 

A 

B 
C 
D 
E 



a s A\7 ad C P T/H^ 
I 



v 3 , 



(V ra d - V a d) ~ ^ V^ 

CD /A 
(2/3) GI 
1/Am 

(2/3) (47rr 2 p) 2 a t A/Am 



(A15) 

(A16) 

(A17) 
(A18) 
(A19) 
(A20) 



1G 



we can give the equation for the specific turbulent kinetic energy in its final form that we implemented in 
our code: 

+C 1 - so)--— + \ E 

On On \ On 

(A21) 

The coefficients A-E are functions of the stellar variables V a d, V ra( j, V M , m, r, p, T, fc rac j, -f/p, Cp, S and y>. 
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Fig. 1. — Core overshoot parameter measured in 
fractions of the pressure scale height for KCT with 
at — 0.11 (circles) and a t — 0.30 (triangles) for 
stellar masses between 1.5 M Q and 15.0 M Q . The 
overshoot from detached eclipsing binaries derived 
by Ribas et al. is given for comparison. KCT sup- 
ports a moderate increase in overshoot with in- 
creasing stellar mass. Overshoot of 0.6 as claimed 
by Ribas et al. at the high mass end cannot be 
explained with KCT. 
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Fig. 2. — Models with no overshoot are in accord 
with observational constraints. Models are found 
in either the (late) main-sequence (dotted, dashed 
lines) or the subgiant phase (solid line). Subgiant 
models can only be found for initial masses less 
than 1.47 M Q and are completely excluded when 
tiny amounts of overshoot of «ov > 0.05 are as- 
sumed. The error box is given by the observa- 
tional constraints on luminosity, radius and effec- 
tive temperature. Its hexagonal shape arises from 
making use of all three constraints. 



21 



CN 
CD 

d 



© 00 

J 00 

J o 
o 



CD 

o 



00 

d 




0.26 0.28 



0.3 



0.32 0.34 0.36 0.38 



Log 1Q (R/Rj 



0- 



Fig. 3. — Best fit overshoot parameters for three 
different overshoot descriptions. The best fit val- 
ues indicate only small overshoot but due to the 
large uncertainties overshoot could range from 
ctOM = . . . 1.18. 
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Fig. 4. — Maximally allowed overshoot parameters 
for different overshoot descriptions when adopting 
the restricted input parameter set (weak limits) 
but still allowing for the full 2a observational error 
box. The model using KCT cannot be shifted into 
the error box despite extremely large values for a t ■ 
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Fig. 5. — Maximally allowed overshoot parameters 
for different overshoot descriptions when adopting 
the full input parameter set (strong limits). The 
model using KCT cannot be shifted into the error 
box despite extremely large values for a t . 
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Fig. 6. — Run of non-local convection term of 
equation (2) for KCT (left panel) and the corre- 
sponding turbulent kinetic energy (right panel). 
Increasing the overshoot parameter a t not neces- 
sarily increases the overshoot extent. 
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Fig. 7. — The averaged large spacings are sensitive 
to the stellar radius as expected. They can be used 
to measure the star's radius more precisely. 
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Fig. 8. — Averaged small p-mode spacings versus 
convective core size for the I = mode and models 
with CT2. 
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Fig. 9. — In comparison to the I — 0, the I = 
1 averaged small p-mode spacings show a better 
correlation with convective core size. 
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Fig. 10. — Calculated g-mode periods as a func- 
tion of radial order n from pulsation analysis for 
first three Vs. For given I, the steepness of the best 
fit to the periods gives the period spacings. Most 
properties of the period spacings can be explained 
by asymptotic theory. 
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Fig. 11. — Brunt- Vaisala frequency N as a func- 
tion of stellar mass for a Procyon A model with 
M = 1.505 M® without overshoot (solid curve). 
To show the relative importance of the molecu- 
lar weight gradient, N is also plotted without the 
contribution from V M (dash-dotted curve). 
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Fig. 12. — Comparison of Brunt- Vaisala, fre- 
quency N of a model with M = 1.497M Q and 
KCT overshoot for two different overshoot param- 
eters at = 0.30 (solid curve) and a t — 0.17 (dotted 
curve). Note that the inner cutoff radius is larger 
in the case of more overshoot. 
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Fig. 13. — Correlation between period spacings of 
g-modes and convective core size for lowest order 
I = 1 mode. 
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Fig. 14. — Correlation between period spacings of 
g-modes and convective core size for I = 2 mode. 
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Fig. 15. — Correlation between period spacings of 
g-modes and convective core size for I = 3 mode. 
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Fig. 16. — Same as Figure 13 but for CT2 convec- 
tion models. 
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Fig. 17. — Same as Figure 14 but for CT2 convec- 
tion models. 



36 



3G 



d 



O 







o 

_° CN 

= d 



o 



12 



l=3 

• NO OVERSHOOT 



CT2: 
CT2: 
CT2: 



a=0.20 
a=0.40 
a=0.60 



j&T 



15 





18 



21 



24 



ri| = 3 (minutes) 



Fig. 18. — Same as Figure 15 but for CT2 convec- 
tion models. 
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Fig. 19. — Same as Figure 13 but for CT3 convec- 
tion models. 
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Fig. 20. — Same as Figure 14 but for CT3 convec- 
tion models. 
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Fig. 21. — Same as Figure 15 but for CT3 convec- 
tion models. 



40 







d 



CN 

6 h 



CN 



8 6 



d 



1=2 



* NO OVERSHOOT 

CT2: a=0.20, 2a error box 

* CT2: a=0.20, 1a error box 

* CT2: a=0.40 





O 



28 



30 



32 



n |=2 (minutes) 

Fig. 22. — Model spreads for full parameter range 
(squares) in comparison to models that are con- 
fined to the lcr observational error box (triangles). 
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Fig. 23. — Model spreads for full parameter range 
(squares) in comparison to models with restricted 
mass (triangles). 
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Fig. 24. — Model spreads for full parameter range 
(squares) in comparison to models with fixed hy- 
drogen content (triangles). 
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Fig. 25. — Model spreads for full parameter 
range (squares) in comparison to models with fixed 
metallicity (triangles) . 
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Fig. 26. — Model spreads for full parameter range 
(squares) in comparison to models with fixed mix- 
ing length (triangles). 
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